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SUMMARY 


A unique, cooperative graduate engineering program has been supported 
for fifteen years under NASA Grant NGR 47-003-052. During the 15-year 
period of support, a spectrum of aeronautics-related research was conducted, 
and full-time graduate students obtained financial support and earned 
masters and doctoral degrees while working in the laboratories at NASA 
Langley Research Center, Hampton, Virginia. An abbreviated listing of 
primary program achievements is shown below. 

• A program of joint student advisement was created involving NASA 
engineers and Old Dominion University engineering faculty. 

• The graduate research participants, enrolled in graduate 
engineering degree programs, were awarded stipend and tuition 
support for a two-year period. 

• Total enrollment in the grant program aggregated forty one (41) 
persons, and produced thirty one (31) graduate degrees (26 Masters, 
5 Ph.D's). 

• Approximately 50% of program graduates took up initial employment 
in or remain in an aeronautics-related enterprise. 

• Research results, documented in theses, reports and papers, were 
produced in the following areas: 

1. Aircraft Design, Aerodynamics, Lift/Drag Characteristics 

2. Avionics: Aircraft Guidance and Control 

3. Fluid Mechanics: Inviscid flow. Boundary Layer and Gas 

Dynamic Studies 

4. Solid Mechanics: Composites, Shells, Beams Under Static, 

Impact and Transient Loading 

5. Instrumentation and Measurement Techniques 
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Thermophysical Properties Experiments 

7. Large Space Structures: Dynamics, Control, Reliability 

Analysis 

8. Earth Orbital Dynamics 

9. Environmental Engineering: Water Qual ity Model ing and 

Atmospheric Radiation Processes 

Of the total grant funds received, 50% went directly to student 
stipends; 12% was for indirect costs of grant management; and, the 
average cost of a graduate engineering degree was $27,000. 

Old Dominion University gained an enhancement of its graduate 
engineering programs during their formative years. 
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INTRODUCTION 

In the early nineteen seventies there was an evident and increasing 
need for joint participation in aeronautical education between universities, 
government, and industry. In 1974 this matter of concern was addressed by 
all three proponents at a joint NASA/University Conference on Aeronautics. 
Various steps were taken to enhance the production of engineering graduates 
who by maintaining great interest and sufficient preparation would assume 
responsible positions in a diverse aeronautics industry. One such coopera- 
tive graduate engineering program operated between the Louisiana State Uni- 
versity in Baton Rouge and the NASA Johnson Space Center. The objective of 
that program was to provide a "hands-on" graduate engineering educational 
experience and to couple faculty, students and NASA engineers. The Old 
Dominion University School of Engineering proposed a similar yet distinct 
program which operated initially from the Virginia Associated Research Cam- 
pus (VARC) in Newport News, Virginia. 

The Aeronautics Graduate Research Program commenced in 1971, with the 
primary goal of engaging students who qualified for regular admission to the 
Graduate School of Engineering at Old Dominion University in a graduate 
engineering research and study program in collaboration with NASA Langley 
Research Center, Hampton, Virginia. Research participants were enrolled in 
the program for two calendar years. During this period, students were pro- 
vided a stipend and tuition expenses under a research and study grant from 
NASA Langley Research Center. Students with diverse baccalaureate prepara- 
tion in both mathematics and engineering were attracted to the program with 
the intent of immersing themselves in a study effort and project work at 
Langley Research Center in aeronautics-related subject matter. The princi- 
pal investigator had the responsibility to solicit nationally for program 


candidates possessing qualifications to enter the graduate engineering de- 
partments. Program participants with aeronautics interests were usually 
selected on the Master's degree level, although several participants, who 
were highly recommended, entered the program at the Ph.D. degree level. 

The majority of the Master's students had little difficulty completing all 
of the requirements, including a thesis, within the two-year tenure. Under- 
standing the program's limited time frame, Ph.D. students and their advis- 
ors usually made further arrangements for continuation and degree comple- 
tion. 

Format and Purposes 

The first year was essentially used for full-time study, although the 
principal investigator sought placement of a student in an appropriate 
branch at NASA Langley Research Center as soon as the student began the 
program. Participating NASA branch Engineers suggested several projects of 
interest as students made initial visits, while the student and his guidance 
committee then decided upon a topic of mutual interest. The second year was 
designated for full-time research except for an occasional course which may 
have been taken to support thesis research or otherwise satisfy degree re- 
quirements. During the two summers that were a part of the time committ- 
ment, the participants were engaged in full-time research at NASA Langley 
Research Center. Graduate participants satisfied course requirements by 
attending lectures at the VARC campus, Newport News, and at the main campus 
in Norfolk; they also met the thesis requirement of the Old Dominion Uni- 
versity graduate engineering programs. 

The primary purpose of the Aeronautics Graduate Research program was to 
engage qualified graduate students in a professional research/design experi- 
ence while simultaneously providing the academic environment for full time 
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study. In the early days, the cooperative program was administered from the 
graduate campus in Newport News, Virginia. Research participants found 
study space, library resources, faculty advisement, and most graduate in- 
struction at the Virginia Associated Research Campus (VARC). At the same 
time, participants found strong research interaction among themselves and 
faculty advisors and NASA research engineers. Objectives of the program, 
defined to serve the interests of selected graduate students, the graduate 
engineering departments at Old Dominion University, and the aeronautics 
research interests at NASA Langley Research Center, were successfully met in 
a number of ways. 

How Objectives Were Met 

In terms of program administration and the technical assistance which 
had been provided, it was most gratifying that NASA Langley Research Center 
engineers and their first line management had found the cooperative program 
appropriate for and consistent with their needs. A good measure of that 
level of satisfaction was the ease with which the principal investigator 
could identify several times more projects and job slots at NASA Langley 
Research Center that needed research assistants, than there were enrolled 
program participants. Top management at NASA Langley Research Center took 
satisfaction in knowing that a supported program, although small in people 
and dollars, produced qualified engineers for selected aeronautics special- 
ties. NASA gained needed technical assistance for projects, and the cul- 
mination was an injection of additional expertise into the aeronautics in- 
dustry. The Old Dominion University School of Engineering was pleased with 
the program because it complemented the basic mission of graduate engineer- 
ing education, combining students, faculty, and practitioners of engineering 
in a productive, learning environment. The opportunity attracted worthy 
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students to a program which then assisted in broadening the base of faculty 
research. Lastly, from the student point of view, a graduate and profes- 
sional level experience was gained. The full-time program participants took 
leave of the closed society of graduate student peers in order to be engaged 
in an actual engineering research and design environment. They traded peers 
in the traditional sense for a realistic communications experience involving 
academic engineers (the faculty), NASA engineers and technicians, and the 
administrators of the program. Of course, their expenses were payed as they 
were being encouraged to complete a graduate engineering degree. Graduate 
research participants, who completed degree requirements, were in an excel- 
lent job-seeking position. 

A joint venture in graduate engineering education is described in this 
report in terms of the successes and problems which were experienced in a 
cooperative program of significant duration. The following sections of this 
final report will cover areas such as student selection and supervision, 
statistical information concerning the graduate participants, research 
accomplishments, and the financial arrangement of the program. 
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GRADUATE PARTICIPANTS AND RESEARCH SUPERVISION 

This chapter provides pertinent data about the origins and progress of 
the research participants. The first section covers student selection and 
advisement followed by the statistical data. 

Student Selection and Advisement 

Students were solicited on a nationwide basis with the main require- 
ments being that candidates had met graduate school admission reqirements, 
had a B.S. in Science or Engineering, and had indicated an interest in grad- 
uate study with an aeronautics emphasis. Candidates for the program must 
have displayed not only acceptable academic credentials, but also behavioral 
traits that indicated a potential for interacting productively with a va- 
riety of advisors and peers, and must have showed poise and self- starting 
capabilities. Solicitation for applicants was conducted each year in the 
winter for candidates who were prepared to enter the two-year program in the 
Fall. By April, awards were made and two-year committments were promised. 
Awardees were selected by the principal investigator, originally at the 
Virginia Associated Research Campus (VARC) , in consultation with the appro- 
priate departmental graduate program directors. Aeronautics research stu- 
dents were then enrolled in a department consistent with their background 
and interest. After being immediately introduced to engineers in several 
branches at NASA/Langley Research Center through contacts made by graduate 
faculty and the program principal investigator, each new participant began 
spending a limited amount of time in the NASA laboratories under the super- 
vision of a NASA engineer, who became a co-chairman of the student's gui- 
dance committee. Concurrently, a faculty member was found who would take an 
interest in the student's topical area for research and study, and was 
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designated as guidance committee co-chairman. The guidance and evaluation 
committee was formally appointed by the participant's departmental chairman 
and had at least three members, a majority of whom were Old Dominion Univer- 
sity faculty members. Faculty advisors became more or less involved with 
the student research project, but the faculty committee co-chairman was 
chosen in order to bring some expertise to the work of the NASA branch where 
the research participant was engaged. The faculty guidance committee co- 
chairman was responsible to the academic department chairman to determine 
that the research participants completed all degree requirements. The NASA 
Advisor was very instrumental in the research/design project planning, and 
was also a full participant in the evaluation of student progress. When 
appropriate, NASA advisors were offered research appointments in the School 
of Engineering. Table 1 is included in order to recognize each graduate 
student's faculty advisor with his asociated university engineering depart- 
ment, and each participating NASA advisor and his associated branch. The 
names of program participants are listed chronologically (see also Table 2). 
Joint advisement was an important element of the program enabling successful 
completion of research tasks and degree requirements. 

Essential features in selecting candidates were the "quality" of the 
applicants and the competition over the years for the relatively few tal- 
ented graduate students interested in aeronautics. Competition was keen for 
requisite applicants; it seemed to be increasingly difficult to lure bacca- 
laureate graduates into full-time graduate study, even with a competitive 
stipend. It was the main problen associated with the administration of the 
grant program, the securing of highly qualified graduate students necessary 
to make the program a success academically. Also, it was important to find 
technically able people for the NASA branches who desired the aid of 
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qualified research assistants. 

A promotional mailing was conducted annually and was directed to chair 
persons of engineering departments around the country in the disciplinary 
areas of mechanical, electrical, civil, aeronautical engineering and en- 
gineering mechanics. The purpose was to attract the most able applicants to 
the graduate engineering research participation program in aeronautics. The 
program administrator found in conversations with other aero-related program 
administrators around the country that the pool of qualified students who 
would enter full-time graduate study with an aeronautics emphasis was small, 
relative to the number of graduate programs competing for these persons. 

The recruiting results, treated in more detail in a later section, demon- 
strates the difficulty experienced in securing student participants best 
suited for the program, a problem likely to continue in the competitive 
market for young engineering talent. 

Program Statistics 

During the 15 year tenure of the Aeronautics Research Participation 
Program, a total of 41 graduate students were enrolled at one time or an- 
other. Since students entered the program on different calendar dates and 
finished their research assignments at their own pace, coordinating with 
their faculty and NASA advisors, the number of active program participants 
changed throughout each academic year. The aggregate of program partici- 
pants each spring semester is illustrated in Figure 1. The graph also, 
shows the number of new awardees enrolled each academic year. In order to 
maintain a "steady state" count in the program, new participants were 
accepted as needed and as funding was available. After an annual recruiting 
campaign (starting in 1974-75) by the principal investigator, offers were 
made and positions were accepted by promising candidates. Figure 2, 
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Figure 1. Program Enrollment Data 







displays the total lumber of applicants, the total number of positions 
offered, and the total number of acceptances for each recruiting period. As 
seen on the graph (of Figure 2), the academic year 1974-75 seemed to be the 
peak year for applicants and offered positions, being 20 and 7 respectively. 
A total of six accepting awardees entered during the academic year 1976-77, 
the largest contingent of research participants to enter at one time. 

Of the 41 program participants, 3 were women and 38 were men. All 
three women, Sophia Ashley, Joanne Walsh, and Joan Gosink, received a Mas- 
ters of Engineering in the two-year time period of the program. The program 
administration was pleased to hear that Ms. Sophia Ashley received a 1985 
Federal Engineer of the Year Award. Ms. Ashley was of Greek origin but 
maintained a USA citizenship. There were seven foreign nationals who parti- 
cipated in the graduate aeronautics research program: five Indians, one 

Nigerian, and one Swede. The remaining 33 participants were of USA citizen- 
ship, although several of those had foreign origins. The program director 
attempted to recruit such candidates as women, minority group USA citizens; 
and, secondarily, foreign nationals for the cooperative aeronautics resarch 
participation program. 

The contrasting backgrounds of the program participants contributed to 
a wide range of academic achievement among the graduates. The participants' 
undergraduate degrees ranged from Bachelors of Science in mechanical, 
electrical, aeronautical, and nuclear engineering to Bachelors of Science in 
physics, mathematics, and geology. Only sixteen out of the 41 overall 
enrollees earned undergraduate degrees from Old Dominion University. In 

part measure this was the result of an effort by the program director to 
bring in qualified graduate students from outside eastern Virginia. 

Principal products of the program were the student graduates who 
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completed their advanced degrees working between the University and 
NASA/Langley Research Center, and then moved on into engineering practice. 
During the active years of the program, thirty-one program participants 
completed their graduate degrees; 26 of these were at the masters level and 
5 at the doctoral level. These persons and their thesis titles are listed 
in Table 2. This table is an accumulative record of all graduates of the 
aeronautics research participation program. Items also listed in Table 2 
are academic accomplishments, academic year of program completion, and the 
specific department and the type of degree earned. As a visual reference. 
Figure 3 depicts the number of program graduates for each academic year. 
Happening in both 1974-75 and 1978-79, the largest contingent of graduates 
to finish in the same academic year was five. In contrast the years 1975-76 
and 1981-82 show no production of graduates from the aeronautics program. 

Dan Stead, a former participant in the aeronautics program, is expected 
to complete his thesis and graduate during the academic year 1986-87 at ODU. 
Mr. Stead was unable to finish degree requirements during the two-year 
program tenure, but is working successfully to complete research in 
experimental fluid mechanics. The remaining 9 of the total 41 participants 
withdrew from the program for either personal reasons or because they failed 
to meet program requirements. The number of awardees who failed to meet 
standard requirements other than personal reasons was only 3. These 
graduate participants were asked to resign after sufficient warning that 
their program success was in danger. The final 6 withdrew for various 
personal reasons, usually early in the program. The principal investigator 
tried to fill those participant vacancies by the beginning of the next 
semester. In summary 31 people graduated with full or partial support from 
the grant program, 1 student is continuing with graduation in sight, and 9 
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*ME - "Master of Engineering" in respective department. 



Table 2. Proaram Graduates under NASA grant NGR 47-003-052 
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Figure 3. Number of Graduates Per Academic Year. 


participants dropped out of the program -- Totaling 41 graduate 
participants. 

Insofar as data is complete at this time. Table 3 displays the initial 
occupational status of thirteen of the program graduates. Many of these 
jobs were entered directly after graduation, therefore they may be 
inaccurate for the present. One object. of the grant program was to 
contribute graduates to the national pool of engineers qualified and 
interested in aeronautics-related practice. As is seen from Table 3, all 
but 4 of the 13 graduates are engaged in some form of aeronautics activity. 
To the extent shown, the Aeronautics Research Participation Program has 
gained a rung of its original objective. It is also interesting to note 
that 4 of the graduates listed in Table 3 gained employment at NASA Langley 
Research Center. 
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Table 3. Occupational Status of Program Graduates 
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RESEARCH ACCOMPLISHMENTS 

Grant funds were allocated for Aeronautics research. The principal 
investigator sought to achieve this objective, but in doing so found that 
the final research product would be determined by (1) the needs of 
companion NASA/LaRC Branches, (2) the already developed directions of 
faculty research, and (3) the capabilities and interests of the graduate 
research participants. The direct research results are the thirty (30) 
theses and several papers, published and presented, which were produced. 
However, the graduate research assistants accomplished in many cases 
additional and peripheral tasks for the NASA hosts, as well as assisting 
faculty advisers towa/d expanding their research programs. At least the 
direct research results can be recorded here. 

While the objective was to produce aeronautical research results, while 
educating graduate students, and the majority of the theses fall in this 
area, other topics were occasionally treated, based on an appropriate match 
of interests between students and their advisers. The three broad research 
areas were aeronautics, aerospace, and environmental engineering. Specific 
research results were achieved in the following categories (number of theses 
in parenthesis): 

1. Aircraft Design, Aerodynamics, Lift/Drag Characteristics 
(3) 

2. Avionics: Aircraft Guidance and Control (3) 

3. Fluid Mechanics: Inviscid flow. Boundary Layer and 

Gas Dynamic Studies (6) 

4. Solid Mechanics: Composites, Shells, Beams Under Static, 

Impact and Transient Loading (5) 

5. Instrumentation and Measurement Techniques (4) 


22 


6. Thermophysical Properties Experiments (12) 

7. Large Space Structures: Dynamics, Control, Reliability 

Analysis (2) 

8. Earth Orbital Dynamics (2) 

9. Environmental Engineering: Water Quality Modeling and 

Atmospheric Radiation Processes (3) 

The aeronautics related categories extend through number 6. above, 
incorporating 23 out of the 30 graduate theses which were produced. A 
statement of the objectives, methods and main results for each of the 
research tasks is contained in the theses abstracts which are included as 
Appendix A of this report. The theses are referenceable documents housed in 
the Old Dominion University Library. 

A number of technical papers were published and presented as a result 
of the work accomplished under the grant. Those papers known to the 
principal investigator, and specifically involving the research 
participants, are listed as Appendix B. The research participants' name is 
underlined in each citation. There were probably additional documented 
works which did not come to the principal investigator's attention. 


23 


BUDGETARY SUMMARY 


Funds were expended over a period of fifteen years to accomplish grant 
program objectives. In fact the original grant award carries a date of 
SeptemDer 1, 1971. It is well to note this was a "non-specific" research 
grant in contrast to the traditional unsolicited proposals and grant mechan- 
ism for research awards. Of course, in the case of the subject grant, there 
was a strong graduate education component coupled with the research objec- 
tives; and, a larger than usual fraction of total grant funds was earmarked 
for graduate research participant stipends. 

Stipends paid to research participants amounted to 63% of total, aggre- 
gate direct cost or 56% of total direct plus indirect costs. These aggre- 
gate amounts are. 

Total direct costs $742,390 

Total indirect costs 96,532 

Total expenditure $838,922 

Another way to view the cost of the program is in terms of the nunber 
of research participants who graduated, the degrees produced. Since Ph.D. 
students and Masters students all were supported for a nominal two-year 
tenure in the program, there is no differentiation of costs on the basis of 
degree level. For thirty-one (31) program graduates, the costs per advanced 
engineering degree are. 

Cost per degree, direct $24,000 

Cost per degree, direct plus indirect $27,000 
The above costs were rounded to the nearest thousand dollars. Certainly on 
the basis of current dollars, this is graduate engineering education at 
bargain prices*. 

tit was occasionally necessary to support by other means some students 
beyond their two-year tenure in the program. Those costs are rolled into 
in the above estimates. 
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Finally, to document in more detail how the grant funds were expended. 
Table 4 shows current, aggregated budget status for the main categories of 
expenditures. 


Table 4. Record of Budgeted Expenditures 



Description 

Amount 

1 . 

Salaries and Wages 

Principal Investigator 
Graduate Students 
Faculty/Research Associates 

143.053.00 

467.218.00 
27,616.00 


Subtotal 

637,887.00 

2. 

Fringe Benefits 

47,411.00 

3. 

Supplies (office/lab) 

9,440.00 

4. 

General Expenses 



Tuition 

Publ ications/subscriptions/report preparation 
Computer charges 

28,585.00 

2.925.00 

1.553.00 


Subtotal 

33,063.00 

5. 

Travel 

14,194.00 

6. 

Equi pment 

395.00 


Total Direct Costs 

742,390.00 


Total Indirect Costs 

96,532.00 


Total Expenditure 

$838,922.00 
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APPENDIX A 


Thesis Abstracts Arranged by the Ni 
Subject Area Categories 


1. Aircraft Design, 
Lift and Drag 
Characteristics 


ABSTRACT 


EXPERIMENTAL STUDY OF DELTA WING LEADING-EDGE DEVICES 
FOR DRAG REDUCTION AT HIGH LIFT 

Thomas Dwight Johnson, Jr. 

Old Dominion University, 1980 
ODU Director: A. Sidney Roberts, Jr. 

NASA Director: Dr. Jim Campbell 

Future fighter aircraft requirements specify efficient supersonic 
cruise and high-g maneuverability at high lift. The slender delta wing 
meets the first requirement but has large lift induced drag increments at 
high lift. One method to alleviate the drag is to control the flow at the 
wing leading edge (LE) by means of small LE devices, so as to maintain loc- 
ally attached flow to higher angles of attack and thus increase the level of 
aerodynamic thrust. 

The devices selected for evaluation were the fence, slot, pylon-type 
vortex generator (VG), and sharp leading-edge extension (SLEE). These de- 
vices were tested on a 60° flatplate delta (with blunt LE) in the Langley 
Research Center (NASA) 7- by 10-foot high-speed tunnel at low-speed and to 
angles of attack of 28*. Balance and static pressure mesuranents were 
taken. 

The results indicate that all the devices had significant drag 
reduction capability and improved longitudinal stability while a slight loss 
of lift and increased cruise drag occurred. 
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ABSTRACT 


ANALYSIS OF SPREAD MULTI-JET 
VTOL AIRCRAFT IN HOVER 

Leroy F. Albang 
Old Dominion University, 1980 
ODU Director: Dr. Tiwari 

NASA Director: Richard J. Margason 

An investigation of vertical takeoff and landing (VTOL) aircraft lift 
losses in hover has been conducted to evaluate a method for a simplified 
test technique. 

Three flat plate models were tested to determine their usefulness in 
predicting hover characteristics by comparing results between plate and 
three-dimensional models. Attempts to correlate the out of ground effect 
results for the plate models with the empirical expression for calculating 
jet induced loads were successful for the models which used engine simu- 
lators producing efflux characteristics similar to the three-dimensional 
model engines. Data for the model using engine simulators with characteris- 
tics unlike those of the three-dimensional model simulators could not be 
correlated in this manner. 

In ground effect, comparisons of induced lift loads indicated correla- 
tions between the plate and three-dimensional models were within 2 percent 

of thrust in the height range dominated by the fountain effect. However, 
small outward deflection of the lift engine exhaust was found to cause a 

decrease in beneficial fountain on the order of 5 to 10 percent of thrust. 
Dashpots proved useful in eliminating data scatter caused by flow-induced 
model vibration. The engine arrangement of a midspan-pod VTOL lift-fan 
transport model showed a strong adverse effect on lift loss due to engine 
deflection. Data obtained for the plate models could be correlated to 
three-dimensional results by the application of a geometrical equivalent 
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height correction factor aH D . The correlation of plate and tunnel models 

e 

indicated that lift losses in ground effect were essentially independent of 
the efflux characteristics for the engine simulators. 
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ABSTRACT 


COMPUTER-AIDED DESIGN OF LIGHT AIRCRAFT 
TO MEET CERTAIN AERODYNAMIC AND STRUCTURAL REQUIREMENTS 

Joanne Lynn Walsh 
Old Dominion University 
ODU Director: Dr. John H. Heinbockel 

NASA Director: W. Jefferson Stroud 

Nonlinear mathematical programming techniques are applied to a particu- 
lar aircraft design problem, that of determining the maximum range design 
for light aircraft subject to certain aerodynamic and structural require- 
ments. Search techniques available in an existing nonlinear optimization 
computer program, AESOP, are applied to the design problem. Seven design 
variables are used to define the wing dimensions, the wing skin thickness, 
the cruise velocity, the amount of fuel carried, and the engine horsepower. 
The influence of different design requirements on the range is observed and 
discussed. 
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2. Avionics: 
Guidance and Control 



ABSTRACT 


STABILITY AND CONTROL OF A SLING LOAD 
USING ACTIVE JETS 

Mark D. Pardue 

Old Dominion University, 1977 
ODU Director: Dr. John W. Stoughton 

NASA Director: Dr. Shaugnnessy 

Helicopter sling load work has been increasing over the past decade, 
and accidents have resulted due to loss of control of the helicopter or 
collision with objects. It is desired to dampen swinging of the sling load 
and thus reduce pilot workload and pilot error. This would have a positive 
effect on accident reduction. 

An active jets control system was designed using velocity feedback. 
Values of feedback gain are presented in graphs for various loads, 
helicopters, and single-cable suspension systems. Damping ratios of 0.7 for 
a wide range of parameters are verified by implementation on a helicopter 
sling load simulation. 
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ABSTRACT 


LOW SPEED HANDLING CHARACTERISTICS OF 
A SUPERSONIC TRANSPORT BY EIGENSTRUCTURE ASSIGNMENT 

Michael Peter Lazar 
Old Dominion University, 1981 
ODU Director: A. Nadkarni 

NASA Director: Walt Oehman 

A synthesis procedure, based upon optimal regulator theory and eigen- 
value/eigenvector assignment technique, was used in the design of a low- 
speed lateral control system for the SCAT-15F supersonic transport. Empha- 
sis was placed on reducing the lateral acceleration at the pilot’s station. 
Control laws were designed to make the aircraft handling characteristics 
acceptable with the exception of the lateral acceleration. It was found 
that a compromise between the roll performance requirements and the peak 
lateral acceleration requirement is needed for a successful SCAT-15F 
design. 
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ABSTRACT 


A COMPARISON OF DIGITAL COMPUTERS IN 
IMPLEMENTING AVIONIC CONTROL LAWS 

Charles C. Marshall 
Old Dominion University, 1977 
ODU Director: Dr. Mural idhara Varanasi 

NASA Director: David R. Downing 

Direct digital control laws impose definite processing requirements 
upon the digital computer. A computer's ability to perform these require- 
ments at a minimum life-cycle cost should be evaluated prior to its 
acquisition. The possibility of using several microprocessors in a multi- 
computer organization to replace a single minicomputer is attractive only if 
it is cost effective; thus the actual life-cycle cost incurred in each 
instance must also oe evaluated. Because of the many interacting effects of 
the possible configurations of a computer systen a straightforward 
evaluation is not feasible. This thesis develops an evaluative tool which 
accounts for these interacting effects. The CPU execution time of one 
iteration of the control law and the affected life cycle costs are combined 
into a figure of merit which can be used to evaluate the relative ability 
and costs of candidate computers. Results indicate that a multicomputer 
organization of microprocessors can perform as well as a single minicomputer 
at one-third the cost. 
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3. Fluid Mechanics: Inviscid, 

Boundary Layer and 
Gas Dynamics Studies 
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ABSTRACT 


AN EXPERIMENTAL INVESTIGATION OF A TURBULENT BOUNDARY LAYER 
WITH SUCTION THROUGH CLOSELY SPACED STREAMWISE SLOTS 

Stephen P. Wilkinson 
Old Dominion University, 1978 
ODU Director: Dr. Robert Ash 

NASA Director: Leonard Weinstein 

The object of this research was to experimentally evaluate a suction 
surface for use with incompressible turbulent boundary layers. The surface 
consisted of an array of closely spaced slots aligned in the direction of 
the free-stream flow. Direct drag and mean boundary- 1 ayer velocity profile 
measurements showed that the slotted surface had nominally the same suction 
characteristics as a porous surface. 
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ABSTRACT 


APPLICATION OF BOUNDARY LAYER THEORY TO SUCTION THROUGH 
STREAMWISE SLOTS IN WIND TUNNEL WALLS 

Michael H. Berger 
Old Dominion University, 1974 
ODU Director: S. N. Tiwari 

NASA Director: I. E. Beckwith 

A theoretical investigation at Mach 6 and 8 has been conducted on a 
model of a wind tunnel sound shield. The flat plate model consists of 
interchangeable circular and wedge shaped longitudinal rods with gaps 
between the rods for boundary layer removal. One-dimensional gas dynamics 
results have been applied to model the pressure distribution for the array 
of cylinders. Comparisons with experimental pressure distributions and 
numerical potential flow solutions indicate good agreement. It has been 
demonstrated that stagnation line heat transfer can be accurately predicted 
for laminar and turbulent flow for circular cylinders. The flow around the 
wedge rods was apparently disturbed and hence no truly laminar or turbulent 
data were obtained. In view of the large effective sweep angle of the 
circular rods with respect to the local flow (89*), boundary layer behavior 
was surprisingly well represented by swept cylinder theory. A study of 
transition data indicates laminar flow was maintained up to a local diameter 
Reynolds number of 1.5 x 10 5 which is in agreement with results for isolated 
swept cylinders. Measurements by NASA - Langley staff indicate considerable 
noise reductions in the shielded region. Thus the model is apparently a 
viable sound shield. 
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ABSTRACT 


EVALUATION OF HYDROGEN AS A CRYOGENIC 
WIND TUNNEL TEST GAS 

Richard Carl Haut II 
Old Dominion University, 1977 
ODU Director: Or. John Kuhlman 

NASA Director: Jerry Adcock 

A theoretical analysis of the properties of hydrogen has been made to 
determine the suitability of hydrogen as a cryogenic wind tunnel test gas. 

By using cryogenic hydrogen, instead of air or cryogenic nitrogen, as the 
wind tunnel test gas, a significant increase in test Reynolds number may be 
achieved without increasing the aerodynamic loads. Under sonic conditions, 
for example, compared to air at ambient temperature, cryogenic hydrogen at a 
pressure of one atmosphere produces an increase in Reynolds number of a 
factor of approximately 14 is while cryogenic nitrogen, at the same 
pressure, produces an increase of only a factor of about 6. 

The theoretical saturation boundary for parahydrogen is well defined. 
Thus, any possible effects caused by the liquefaction of the test gas can 
easily be avoided by knowing the maximum local Mach number on the model. 

The nondimensional ratios used to describe various flow situations in 
hydrogen were determined and compared with the corresponding ideal diatomic 
gas ratios. The results were used to examine different inviscid flow con- 
figurations. This investigation concluded that the relatively high value of 
the characteristic rotational temperature causes the behavior of hydrogen, 
under cryogenic conditions, to deviate substantially from the behavior of an 
ideal diatomic gas in the compressible flow regime. Therefore, if an ideal 
diatomic gas is to be modeled, cryogenic hydrogen is unacceptable as a wind 
tunnel test gas in a compressible flow situation. 

However, at low Mach numbers where the assumption of incompressiblity 


39 


is valid, the deviations in the isentropic flow parameters for cryogenic 
parahydrogen from the corresponding flow parameters for an ideal diatomic 
gas are negligible. Thus, in the incompressible flow regime, cryogenic 
hydrogen is an acceptable test gas. 

Hydrogen properties and fan drive-power requirements related to a hy- 
drogen wind tunnel were also examined. The drive-power requirements were 
found to decrease wth decreasing temperature and may be adequately predicted 
by using modified versions of the ideal gas equations. Since gaseous hydro- 
gen is capable of penetrating and degrading the mechanical characteristics 
of numerous materials, materials known to be compatible with hydrogen must 
be used exclusively in the design of a hydrogen wind tunnel to avoid prob- 
lems as a result of exposure to gaseous hydrogen. 

A literature survey resulted in the conclusion that although hydrogen 
is a highly combustible substance, safety codes exist which, when followed, 
minimize the risk involved in handling hydrogen. 
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ABSTRACT 


INVISCID TWO-DIMENSIONAL VORTEX-PANEL METHOD FOR CALCULATING 
THE PRESSURE DISTRIBUTION OVER NONCIRCULAR CYLINDERS 
AT VARIOUS FLOW INCIDENCE ANGLES 

Sophia Kokkinos Ashley 
Old Dominion University, 1978 
ODU Director: Dr. John M. Kuhlman 

NASA Director: Ed Polhamus 

A method has been implemented to calculate the potential flow pressure 
distribution over noncircular cylinders at various flow incidence angles. 

The inviscid solution has been obtained by replacing the surface of a 
noncircular cylinder with a vortex sheet. The surface has been approximated 
by a number of flat panels. Each vortex panel has been assumed to have a 
constant vorticity, and the condition of no-circulation applied. Solutions 
yield the local tangential velocity and pressure coefficient at each point. 
The accuracy of the method has been established through comparison of the 
predicted results with existing exact theoretical results. 

Comparisons of the predicted inviscid results with experimental data 
for square and triangular cross-sectional cylinders with rounded corners 
have also been included. Details of the computer program developed to 
implement this method are described in Appendixes A and B, which include 
input and output examples for a sample problem. 
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ABSTRACT 


INTERACTION OF A TWO-DIMENSIONAL TRANSVERSE 
JET WITH A SUPERSONIC MAINSTREAM 

Gilbert 0. Kraemer 
Old Dominion University, 1982 
ODU Director: Dr. Surendra N. Tiwari 

NASA Director: Dr. R. Clayton Rogers 

An experimental investigation of the interaction of a two-dimensional 
sonic jet injected transversely into a confined main flow has been con- 
ducted. The main flow consisted of a air at a Mach number of 2.9. The 
effects of varying the jet parameters on the flow field were examined using 
surface pressure and composition data. Also, the downstream flow field was 
examined using static pressure, pitot pressure, and composition profile 
data. The jet parameters varied were gapwidth, jet static pressure, and 
injectant species of either heliun or nitrogen. The values of the jet pa- 
ramaeters used were 0.039, 0.056, and 0.109 cm for the gapwidth and 5, 10, 
and 20 for the jet to mainstream static pressure ratios. The features of 
the flow field produced by the mixing and interaction of the jet with the 
mainstream were related to the jet momentum. The data were used to demon- 
strate the validity of an existing two-dimensional elliptic flow code. 
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ABSTRACT 


NUMERICAL STUDY 

OF THE TURBULENT BOUNDARY LAYER 
IN A HYPERSONIC NOZZLE 

Joan Patricia Gosink 
Old Dominion University 
ODU Director: Dr. G. L. Goglia 

NASA Director: Julius E. Harris 

A numerical study of the turbulent boundary layer in a hypersonic noz- 
zle is presented. The boundary layer equations are solved by a second or- 
der, implicit finite difference scheme developed at Langley Research Center. 
The technique is applicable to laminar, transitional or turbulent compres- 
sible or incompressible boundary layer flows for planar or axisymmetric 
geometries. A technique is presented for utilizing experimental data to 
start the finite difference solution. 

Turbulent boundary-! ayer measurements on axi symmetric nozzle walls at 
Mach Number 6 are presented and compared with numerical predictions obtained 
by starting the solution in the stagnation chamber and from a downstream 
experimental station. Accurate numerical predictions were obtained for all 
cases when the classical boundary layer equations were valid. The numerical 
solution utilizing local data to start the finite difference solution 
improves agreement with measured velocity improves agreement with measured 
velocity and temperature profiles and displacement thickness. 

Boundary layer measurements and numerical results indicate that the 
Crocco-velocity relation relaxes from a quadratic- type variation toward a 
more linear variation as the pressure gradient relaxes to zero. 
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For low Reynolds number flows modification of the eddy viscosity model 
utilized in the present finite difference scheme improves the accuracy of 
the predictions. Local application of an approximate low Reynolds mmber 
eddy viscosity model indicates improved agreement in boundary layer 
thickness, displacement thickness, and momentum thickness. 
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4 . 


Solid Mechanics: 
Static, Impact, 


Composites, Shells, Beams, 
and Transient Loading 


ABSTRACT 


AN ANALYTICAL ANO EXPERIMENTAL STUDY OF CLEARANCE AND BEARING- 
BYPASS LOAD EFFECTS IN COMPOSITE BOLTED JOINTS 

Rajiv Vikas Arun Naik 
Old Dominion University, 1986 
Director: Dr. Ram Prabhakaran 

Co-Director: Dr. John H. Crews, Jr. 

A combined analytical and experimental study is conducted to determine 
the effects of clearance and bearing-bypass loading for mechanically 
fastened joints in composites. A simple method of analysis is developed to 
account for the nonlinear effects of bolthole clearance. The nonlinear 
load-contact variations for clearance-fit fasteners are also measured using 
specially instrumented fasteners. For a quasi-isotropic graphite/epoxy 
laminate, results show that the contact arc as well as the peak stresses 
around the hold and their locations are strongly influenced by the 
clearance. After a slight initial nonlinearity, the peak stresses vary 
linearly with applied load. The typical clearance levels are shown to have 
only a minor influence on the overall joint stiffness. 

Quasi-isotropic graphite/epoxy lami nates (T300/5208) are tested under 
combined bearing and bypass loading to study failure modes and strengths. 
Radiographs are made after damage onset and after ultimate load to examine 
the failure modes. Also the laminate stresses near the bolt-hole are 
calculated for each test condition, and then used with appropriate failure 
criteria to analyze the test data. The tension data show a linear 
interaction for combined bearing and bypass loading with damage developing 
in the net-tension mode and growing to failure in the same mode. Failure 
modes are more complex in compression. The compression bearing-bypass 
strangths for damage-onset show an unexpected interaction involving the 
bearing mode. Compressive bypass loads reduce the bearing strength by 
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decreasing the bolt-hole contact arc and thus increasing the severity of the 


bearing loads. Bearing damage-onset, for compressive bearing-bypass loads, 
causes a weakening of the offset-section leading to ultimate failure in the 
offset compression-bearing mode. Damage-onset is predicted reasonably well 
using the peak stresses at the hole. Strength predictions indicate that 
damage corresponding to ultimate strength is governed by the maximum stress 
near the hole. 

This investigation will help improve the basic understanding of 
composite bolted joints and lead to better structural design procedures. 
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ABSTRACT 


ANALYSES OF QUASI- ISOTROPIC COMPOSITE PLATES UNDER QUASI-STATIC 
POINT LOADS SIMULATING LOW-VELOCITY IMPACT PHENOMENA 

Ajit Dhundiraj Kelkar 
Old Dominion University, 1985 
ODU Director: Or. Ram Prabhakaran 

NASA Director: Dr. W. Elber 

Composite laminates have high strength to density ratios that make them 
attractive for use in aircraft structures. However, the damage tolerance of 
these materials is limited because they have very low ultimate strains, no 
plastic deformation range, and no usable strength in the thickness 
direction. These limitations are very obvious when laminates are subjected 
to impact loads. Due to these impact loads, laminates suffer visible and 
invisible damage. To improve the material performance in impact requires a 
better understanding of the deformation and damage mechanics under impact 
type loads. 

In thin composite laminates, the first level of visible damage occurs 
on the back face and is called "back face spalling." A plate-membrane 
coupling model, and finite element model to analyze the large deformation 
behavior of eight-ply quasi-isotropic circular composite plates under impact 
type point loads are developed. The back face spalling phenomenon in thin 
composite plates is explained by using the pi ate- membrane coupling model and 
the finite element model in conjunction with the fracture mechanics 
principles. The experimental results verifying these models are presented. 
The study resulted in the following conclusions: 

1. The large deformation behavior of circular isotropic membranes 
subjected to arbitrary axi symmetric loading can be obtained by 
solving a single nonlinear governing equation in terms of a radial 
stress. 
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2. Accurate large deflection behavior of circular quasi-isotropic 
T300/5208 laminates can be obtained by using a simple plate- 
membrane coupling model. 

3. The functional form of deformed shape of the plate undergoing large 
deformations is different from the small deflection plate 
solution. 

4. The back face spalling action in thin composite laminates is a 
spontaneous action and can be predicted by using the fracture 
mechanics principles. 

5. Mixed mode (I and II) type deformations probably occur during back 
face spalling, however, mode I appears to govern the delamination 
growth during the spalling action. 
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ABSTRACT 


THE EFFECT OF UNREINFORCED CUTOUTS ON THE BUCKLING OF THIN 
CONICAL SHELLS LOADED BY CENTRAL AXIAL COMPRESSION 

Garla Arjuna 

Old Dominion University, 1974 
NASA Director: Dr. James H. Starnes 

ODU Director: Dr. Kenneth H. Murray 

Analytical trends for the linear bifurcation buckling of thin isotropic 
right-circular-conical shells with unreinforced cutouts are investigated. 
Clamped, constant-thickness conical shells with homogeneous, isotropical 
material properties are analyzed for uniform axial compression. The 
buckling loads are obtained by the use of the STAGS STructural Analysis of 
General Shells) computer program which is based on a finite difference 
numerical approach. Variables defining the geometry of the cutout and the 
cone are nondimensional ized and the relationship between the nondimensional 
buckling load and nondimensional cutout parameter is established. A range 
of the cut out sizes is investigated for which linear bifurcation buckling 
results are expected to be reasonable approximations to the actual buckling 
loads. 
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ABSTRACT 


AN EXPERIMENTAL AND ANALYTICAL INVESTIGATION 

OF THE 

IOSIPESCU SHEAR TEST FOR COMPOSITE MATERIALS 

Barry Stuart Spigel 
Old Dominion University, 1984 
ODU Director: Dr. R. Prabhakaran 

NASA Director: Dr. James W. Sawyer 

Mechanical properties of composite materials under shear loading are 
difficult to determine. The Iosipescu Shear Test, originally proposed for 
metals, has in recent years been applied to composites. It has the advan- 
tages of small specimen size, simple loading and a reasonably uniform snear 
stress in the test section. 

The purpose of this work is to study the validity of tne Iosipescu test 
method for measuring the shear modulus and shear strength of composites. 
Finite element analyses indicate that optimum specimen geometry and load 
locations depend upon the degree of orthotropy of the composite. Test re- 
sults for a quasi-isotropic graphite/epoxy laminate show that modifications 
in the notch geometry aimed at improving stress distributions result in 
unexpected changes in the failure mode. Thus, while the Iosipescu Shear 
test is valid for some composites with a suitable choice of test parameters, 
its validity for other composites, especially with a higher orthotropy, is 
doubtful . 


i 
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ABSTRACT 


FINITE ELEMENT MODEL OF A TIMOSHENKO BEAM 
WITH STRUCTURAL DAMPING 

Louis Ablen Roussos 
Old Dominion University, 1980 
Director: Earl A. Thornton 

NASA Director: Dr. E. T. Kruszewski 

A numerical integration technique, a modified version of the Newmark 
method, is applied to transient motion problems of systems with mass, stiff- 
ness, and small nonlinear damping. The nonlinearity is cast as a pseudo- 
force to avoid repeated recalculation and decomposition of the effective 
stiffness matrix; thus, the solution technique is dubbed the "pseudo- force 
Newmark method." Comparisons with exact and perturbation solutions in sin- 
gle-degree-of- freedom problems and with a Gear-method numerical solution in 
a cantilevered Timoshenko beam finite element problem show the solution 
technique to be efficient, accurate, and, thus, feasible provided the non- 
linear damping is small. As a preliminary step into the investigation of 
the active control of large space structures, a problem involving a free- 
free Timoshenko beam with nonlinear structual damping is solved. As ex- 
pected, small damping is shown to be of little importance in the prediction 
of low-frequency vibrations while being of utmost importance in the predic- 
tion of high-frequency vibrations. 
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5. Instrumentation 
and Measurement 
Techniques 


ABSTRACT 


STATIC TEMPERATURE MEASUREMENTS IN SUPERSONIC FLOW 
USING LASER RAMAN SPECTROSCOPY 

Paul A. Nystrom 
Old Dominion University, 1977 
ODU Director: A. S. Roberts- Jr. 

NASA Director: Dr. J. M Seiner 

Laser Raman spectroscopy is a nonintrusive technique which can generate 
mean static temperature (MST) data within flow fields. By ratioinq two 
simultaneously gathered portions of the Raman scattered spectrum, absolute 
intensity measurements were not needed. This method was applied to a 
supersonic cold air iet, M = 1.99, where the MST ranged from 160 # K to 300°K. 
Accurate measurements were not realized due to the presence of flow 
contaminates, especially water vapor. 


ABSTRACT 


A HIGH-SPEED ALGORITHMIC SIGNAL PROCESSOR 

Stephen P. Livingston 
Old Dominion University, 1978 
ODU Director: Dr. Mural idhara Varanasi 

NASA Director: Mr. Patrick Gaines 

Microprocessors are now being used in the implementation of solutions 
to many classes of problems, but in some cases the microprocessor is not 
sufficient due to the lack of processing speed of large-scale integrated 
circuitry. A possible solution is presented using a hardware algorithmic 
signal processor to take the place of time-consuming software routines, 
leaving the microprocessor the task of controlling operations within the 
system. The hardware processor consists of a high-speed multiplier as the 
central unit and high-speed Schottky logic, all under microprogram control. 
Several classes of algorithms are considered for implementation as a 
microprogram for the evaluation of transcendental functions and floating 
point arithmetic using the multiplier as the central unit. Design 
considerations include high speed, which will be evaluated for the systen, 
and small size, which means a decrease in integrated circuit package count. 
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ABSTRACT 


SECOND ORDER OPEN LOOP OMEGA PHASE DETECTOR: 
MICROPROCESSOR-BASED 

Phillip S. Mastic 
Old Dominion University, 1978 
ODU Director: Dr. Preston B. Johnson 

NASA Director: Mr. Carroll Lytle 

A microprocessor-based open loop OMEGA phase detector has been 
investigated and experimentally verified. The microprocessor's ability to 
perform the required phase detection algorithms in software drastically 
reduces hardware complexity providing an OMEGA navigation receiver design of 
lower cost. Also, algebraic manipulation of the quantized phase error by 
the processor linearizes the phase detector characteristic curve. 

This thesis describes the theory and implementation of this open loop 
technique and a comparison is made with a standard digital phase-locked 
loop. A brief review of the OMEGA navigation system is presented followed by 
a discussion of the current receiver developed by personnel at NASA Langley 
Research Center. Succeeding chapters document theoretical and experimental 
efforts made in support of tne goal to upgrade the phase detector section in 
the current NASA receiver. 
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ABSTRACT 


THE DEVELOPMENT OF A MICROPROCESSOR 
CONTROLLED LINEARLY ACTUATED VALVE ASSEMBLY 

Raymond H. Wall 
Old Dominion University, 1984 
ODU Director: Dr. Robert L. Ash 

NASA Director: Dr. Ping Tcheng 

The development of a proportional fluid control valve/assembly is pre- 
sented. This electro-mechanical system is needed for space applications to 
replace the current proportional flow controllers that use bulky hydraulic 
or pneumatic actuation systems. The flow is controlled by a microprocessor 
system that monitors the control parameters of upstream pressure and 
requested volumetric flow rate. The microprocessor achieves the proper 
valve stem displacement by means of a digital linear actuator. A linear 
displacement sensor is used to measure tne valve stan position. This 
distance is monitored by the microprocessor system as a feedback signal to 
make the control closed-loop. With an upstream pressure between 15-47 psig, 
the developed system operates between 779 standard cm 3 /sec (SCCS) and 1543 
SCCS. The delivered volumetric flow rates were within 5% of the requested 
values. A statement on the system stability is included as an appendix. 
Recommendations made for future work on the design include further 
modification to the bellows spring and the possibility of a mass flow rate 
controller. 
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6. Thermophysical Properties 
Experiments 
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ABSTRACT 


A METHOD FOR IMPROVING THE ACCURACY IN PHASE CHANGE 
HEAT TRANSFER DATA THROUGH INCREASED PRECISION IN 
THERMOPHYSICAL PROPERTY DETERMINATIONS 

John Philip Drummond 
Old Dominion University, 1975 
ODU Director: Dr. Robert L. Ash 

NASA Director: Mr. Robert A. Jones 

Wind tunnel model materials typically have thermophysical properties 
that vary with both position and temperature. A numerical study using a 
finite difference representation of the heat conduction equation has been 
performed on several of these model materials to determine their thermal 
response under a specified heat loading. This study has confirmed the fea- 
sibility of using an effective thermophysical property to account for the 
existance of variable thermal properties when a model undergoes aerodynamic 
heating. The effective property, when used in the reduction of phase-change 
heat transfer data, forced agreement of the actual transient behavior of a 
model with the constant property semi-infinite slab approximation tradition- 
ally used to describe model response. Results of the numerical study have 
also been compared with similar data from an experimental apparatus cur- 
rently being developed at NASA Langley Research Center. Knowledge gained 
from the comparison has been used to formulate methods for improving the 
performance of the experimental apparatus. 
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ABSTRACT 


TEMPERATURE AND TIME DEPENDENT VISCOSITY 
OF AN ELASTOMER-MODIFIED EPOXY RESIN 

Woodrow W. Wagner, Jr. 

Old Dominion University, 1978 
ODU Director: Dr. John K. Kuhlman 

NASA Director: Howard Price, Jr. 

Temperature effects on the viscosity of an epoxy resin (diglycidyl 
ether of bisphenol A - DGEBA) , an elastomer (carboxy-terminated butadiene 
acrylonitrile - CTBN), and two blends nominally containing 10 percent elas- 
tomer were examined from strain rates of 0.1 sec- 1 to 10,000 sec- 1 over a 
temperature range from 300K to 450K. All results indicated that the great- 
est drop in viscosity occurred for each of the materials over the tempera- 
ture range of 300K to 350K. In addition, the materials exhibited Newtonian 
responses over the entire strain rate range. Adding small amounts of elas- 
tomer to the epoxy only moderately affected the viscosity; whereas the cata- 
lyst, tri phenyl phosphine, added to one of the blends did not measurably 
change the viscosity. Activation energies over the test temperature range 
were calculated using the Arrhenius relationship between viscosity and in- 
verse temperature. Two distinct activation energies were found for eacn 
material in the temperature ranges 300K to 350K and 400K to 450K, indicating 
the existence of two different molecular flow mechanisms in these two tem- 
perature ranges. 

An elastomer-modified expoxy resin and the epoxy resin by itself were 
tested at constant strain rates as a function of time under conditions of 
room temperature (300K) . Properties of thixotropy and rheopoxy were found, 
observed, and analyzed. Each particular strain rate contained its own char- 
acteristic relaxation and recovery times for a particular material. The 
elastomer showed no time-dependent effects up to a strain rate of 100 sec -1 ; 
reliable data could not be obtained at higher strain rates. 
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Large Space Structures: Dynamics, 

Control, Reliability Analysis 
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ABSTRACT 


OPTIMAL CONTROL OF A LARGE SPACE TELESCOPE 
USING AN ANNULAR MOMENTUM CONTROL DEVICE 

Arun Anant Nadkarni 
Old Dominion University, 1977 
ODU Director: Dr. W. J. Breedlove, Jr. 

NASA Director: Dr. W. W. Anderson and Mr. Nelson J. Groom 

Application of a new development in the field of momentum storage de- 
vices, the Annular Momentum Control Device (AMCD), to the twin problems of 
large angle maneuvers and fine pointing control is considered. The basic 
concept of AMCD consists of a spinning rim, with no central hub area, sus- 
pended by a minimum of three magnetic bearings, and driven by a noncontact- 
ing electromagnetic spin motor. The dissertation considers in detail the 
application of the AMCD to achieve both the large angle maneuvers and the 
fine pointing control of a Large Space Telescope (LST) with a single 
configuration, consisting of a single AMCD Mounted in a single gimbal. 

The problem of designing the optimal controller is accomplished in two 
parts: 1) the optimal controller for generating the open-loop control law 

for the nonlinear maneuvering problan was designed using a modified gradient 
technique with penalty functions, 2) the optimal stochastic controller for 
generating the constant gain, state feedback control law for the linear fine 
pointing problem was designed. 

The open-loop optimal control for the maneuvering problem was derived 
iteratively, using a modified gradient technique with penalty functions. 
"Hard" constraints on gnitudes of state and control variables were incor- 
porated in the gradient procedure. A general, user-oriented computer pro- 
gram was developed to derive open-loop optimal control law history using 
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this procedure for a general high order (up to 25) nonlinear systan. The 
program was used to achieve the design objective of deriving the open-loop 
optimal control law to retarget the LST from one stellar target to another 
with the expenditure of minimum energy. A specific example problem of man- 
euvering the LST through prescribed pointing angles was solved. Convergence 
of the iterative procedure to a local minimum was shown to be highly depen- 
dent on the intitial control history chosen and choice of weighting mat- 
rices. 

The optimal state feedback control law for the fine pointing controller 
was derived using Linear-Quadratic-Gaussian (LQG) optimal regulator theory. 
Structural ly, this controller consists of an optimal regulator and a Kalman- 
Bucy filter. Existing techniques to compute the initial stabilizing gains 
and iterative procedures to solve the algebraic Riccati equation were dis- 
cussed. Extension of these techniques to the present problem, which is 
uncontrollable, was indicated by discussing the modifications of tne state 
equations required to make the regulator problem stabilized. It was shown 
that the initial errors in the states can be nulled satisfactorily. The 
performance of the fine pointing controller was investigated by performing a 
linear covariance analysis to obtain the RMS pointing errors. The analysis 
indicated that fine pointing accuracies of less than one arcsecond could be 
achieved. 
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ABSTRACT 


AN APPROACH FOR INCLUDING RELIABILITY 
IN CONTROL SYSTEM DESIGNS 

Ime Isaac Akpan 
Old Dominion University, 1982 
ODU Director: A. Sidney Roberts, Jr. 

NASA Director: Raynond C. Montgomery 

An approach is developed for non-catastrophic control component failure 
accommodation for large space structures. Reliability is considered for 
both the control configuration (i.e., selecting actuator locations over an 
admissible set) and the ultimate operation of the system. The approach 
presented is to construct a measure of performance that indicates the abso- 
lute goal of the system, and for a given structural design, to determine the 
lowest achievable cost conditional upon the failure states of the system 
(i.e., which actuators are operational). These costs are then weighted and 
summed to provide an overall system performance measure. This expected 
performance measure is then minimized over the set of admissible control 
configurations. A complete mathematical development is presented for a 
flexible beam and a grid structure. The optimal configurations for the 
failure states are determined for different normalized time of mission. 
Although the research is concerned with a particular problem, the analysis 
of the reliability incorporation in the design of control systems is of more 
general interest. 
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8. Earth Orbital Dynamics 
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ABSTRACT 


DYNAMICS OF EXPLOSION REMNANTS IN EARTH-ORBITS 
Joseph T. Fuss 

Old Dominion University, 1974 
NASA Director: Mr. J. M. Alvarez 

ODU Director: Dr. W. J. Breedlove, Jr., 

Explosions in space produce particles that are hazardous to space vehi- 
cles. Some of these particles are not trackable and their locations must be 
predicted by dynamical methods. 

This investigation sheds light on three different areas. First, ignor- 
ing perturbations, the characteristics of a three dimensional envelope out- 
lined by the sum of individual particle orbits is determined. Second, the 
distribution of the particles is found as a function of time. Finally, the 
movement of the body via secular perturbations is investigated. Magnitudes 
of the explosion, the point of detonation in the parent satellite's orbit, 
and the radius of perigee of the parent satellite are varied as parameters. 

It was found that the particles evenly distributed themselves through- 
out the envelope, with areas of higher density due to its shape. Perturba- 
tions caused drastic changes in the size and orientation of the envelope, 
thereby complicating the prediction of the areas of higher density. 
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ABSTRACT 


DYNAMIC ANALYSIS OF AN EARTH-ORBITING SATELLITE 
DEPLOYING HINGED APPENDAGES 

Kevin Fielding Oakes 
Old Dominion University, 1986 
ODU Director: Dr. Chuh Mei 

NASA Director: Dr. J. M. Housner 

The planar librational dynamics of an arbitrary spacecraft with 
deployable hinged appendages is presented here both analytically and 
numerically. The system model consists of a central rigid satellite moving 
about the earth in a circular orbit at a uniform orbital rate. Attached to 
the satellite are two unfolding rigid booms with specified mass and length. 
Tip masses are included at the ends of each boom and deployment is driven by 
torsional springs at both frictionless hinge points. The equations of 
motion are developed for both a single- and double-appendage model through 
the use of Lagrange's equations. The formulation accounts for large angle 
rotations, Coriolis effects, and the gravitational gradient. Utilizing 
Newtonian mechanics, the equations of motion for the deployment of a single 
massless boom are developed and compared with the Lagrangian formulation for 
validation purposes. The resulting nonlinear, coupled governing equations 
are not subject to closed-form solutions ana are thus integrated 
numerically. Nondimensional parameters are defined to simplify the single 
appendage analysis and the equations are linearized about prescribed 
positions to study the stability of the system. The investigation is 
applied to the Space Shuttle based deployment of rigid truss-like members. 
Results indicate that spacecraft inertia parameters, appendage mass and 
length, deployment velocity, and initial conditions all exhibit considerable 
influence on the system response. The resulting librational motion is 
directly related to the size of the deployable payload and gravitational 
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forces lead to vehicle stabilization. These findings may be helpful in 
assessing orbital effects during the deployment of large space structures 
on the motion of satellites with robotic manipulator arms. 
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ABSTRACT 


A MODEL TO INVESTIGATE THE INFLUENCE 
OF SUSPENDED SEDIMENT ON THE MASS TRANSPORT 
OF A POLLUTANT IN OPEN CHANNEL FLOW 

Raymond Scott Chapman 
Old Dominion University, 1977 
ODU Director; Dr. Chin Y. Kuo 
NASA Director: Charles Whitlock 

The environmental impact of the transport of pollutants in open channel 
flow has for many years been of interest due to the continuous introduction 
of heavy metals, pesticides, herbicides, and other foreign substances into 
natural waterways. In order to fully understand this transport process, it 
is necessary to examine the significance of its individual components. In 
the present study an explicit two-dimensional finite difference model, de- 
signed to investigate the influence of suspended sediment on the pollutant 
transport process, is presented. Specific attention is directed toward 
examining the role of suspended sediment in 1) the turbulent vertical trans- 
port mechanism in a stratified flow, and 2) pollutant uptake due to sorp- 
tion. Results presented indicate that suspended sediment plays a major role 
in tne pollutant transport process, and subsequently, any meaningful 
attempt to model the fate of a pollutant in an alluvial channel must account 
for the presence of a suspended sediment concentration profile. Similarly, 
the vertical and longitudinal pollutant concentration distributions provided 
by the model may be utilized to improve upon the predictive capacities of 
existing water quality models. 
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EFFECTS OF MULTIPLE SCATTERING AND SURFACE 
ALBEDO ON THE PHOTOCHEMISTRY OF THE TROPOSPHERE 

Tommy Reinhold Robert Augustsson 
Old Dominion University, 1981 
ODU Director: Dr. Surendra N. Tiwari 

NASA Director: Dr. Joel S. Levine 

A one-dimensional photochemical model of the troposphere containing the 
species of the nitrogen, oxygen, carbon, hydrogen, and sulfur families has 
been developed and used to investigate the vertical profiles and the natural 
(including atmospheric chemical reactions) and anthropogenic sources and 
sinks of these species. The species continuity equations are solved numeri- 
cally applying prescribed boundary conditions. The vertical flux is sim- 
ulated by use of the parameterized eddy diffusion coefficients. Heterogen- 
eous losses (e.g. rainout, gas-to-particle cnemistry, and dry deposition), 
are parameterized to make calculated profiles consistent with measurements. 
The photochemical model is coupled to a radiative transfer model that calcu- 
lates the radiation field due to the incoming solar radiation which initi- 
ates much of the photochemistry of the troposphere. Comparisons of vertical 
profiles of tropospheric species are made between the Leighton approxima- 
tion, widely used in most tropospheric models, and the detailed radiative 
transfer matrix inversion model used in this study. The radiative transfer 
code includes the effects of multiple scattering due to molecules and aero- 
sols, pure absorption and surface albedo on the transfer of incoming solar 
radiation. The results indicate that significant differences exist for 
several key photolysis frequencies and species number density profiles oe- 
tween the Leighton approximation and the profiles generated with the more 
detailed radiative matrix inversion technique used in this study. Most 
species snow enhanced vertical profiles when the more realistic treatment of 


71 


the incoming solar radiation field is included. Furthermore, most species 
increase in concentration as a function of increasing surface albedo. A 
few species, notably ozone, exhibit decreased levels of concentration when 
the realistic radiative transfer model is used. The effect of the detailed 
treatment of incoming solar radiation on the photochemistry of the tropo- 
sphere is discussed. 
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EVALUATION OF UPWELLLING INFRARED RADIANCE 
IN A NONHOMOGENEOUS NONEQUILIBRIUM ATMOSPHERE 

Sundaresa Venkata Subramanian 
Old Dominion University, 1976 
ODU Director: Dr. S. N. Tiwari 

NASA Director: Dr. H. G. Reichle 

The primary purpose of this study was to investigate the influence of 
vibrational nonequilibrium upon upwelling infrared radiance from the Earth's 
atmosphere. Basic equations are presented for calculating the atmospheric 
transmittance, heating rates, and equilibrium and nonequilibrium upwelling 
radiances. By employing the Lorentz line-by-line model for spectral absorp- 
tion, heating rates and upwelling radiances were calculated for equilibrium 
as well as nonequilibrium conditions in the spectral range of 4.6 y CO and 
3.3 y CH 4 bands. For the spectral range of CO fundamental band, the influ- 
ence of different parameters on the upwelling radiance was also investi- 
gated. The results indicate that for CO the assumption of LTE is not justi- 
fied at tropospheric pressures and temperatures. For CHi*, however, the 
assumption is justified up to 60 kilometers. This information is very use- 
ful in developing an accurate data reduction scheme for the measurement of 
CO and CH 4 concentrations in the atmosphere. 
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